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ABSTRACT 

The longitudinal shower development of EAS observed 
in the Fly's Eye is used to determine the distribution of 
X max , the depth in the atmosphere of the EAS maximum. 

Work in progress to compare data and Monte Carlo simulations 
of proton and iron primaries is described. Preliminary 
evidence is in favor of a substantial contribution from 
light primaries. 

1. Introduction. The overall longitudinal development of the EAS 
detected by the Fly's Eye(l) can be used to determine the depth in the 
atmosphere in gm/cm 2 of the EAS maximum X,p ax . The distribution in X max 
is in principle sensitive to the composition of the primary particles 
since iron nuclei and protons will give rise to X max distributions that 
peak at shallower and deeper depths and have narrower and wider widths, 
respectively. It also follows that a mixed composition will have a 
broader X max distribution than any single source. 

In what follows, we discuss the reconstruction of longitudinal 
shower profiles and the systematics of determining X max distributions 
and then discuss work in progress on Monte Carlo simulations, which 
include the details of the Fly's Eye acceptance, for pure protons, 
iron, and a mixed composition. 

2. Shower Size Measurement . A fit to the relative time of arrival 

of light to succeeding phototubes in the event-detector plane yields 
Rp, the impact parameter of the shower to the detector, and the zenith 
and azimuthal angles of the EAS. Measured values of optical gathering 
power, efficiencies and electronic gains and pedestals are used to 
convert photoelectron yields into apparent brightness, i.e., numbers of 
photons arriving at the detector from the source. This can be converted 
into intrinsic source fluorescent brightness after correcting for: (a) 

directly produced Cerenkov light beamed in the direction of the detector; 
(b) Cerenkov light scattered in the direction of the detector due to 
Rayleigh and Mie scattering; and (c) atmospheric attenuation of light. 

The details of these corrections are described in reference 1. The 
intrinsic fluorescence brightness can be translated directly into a 
shower size using the known nitrogen fluorescence efficiency. 

The Cerenkov light production model, and in particular, the 
dependence of the Cerenkov light intensity on emission angle has been 
checked by using a sample of events seen by both Fly's Eye I and Fly's 
Eye II, a smaller station situated 3.3 km from Fly's Eye I. We find 
good consistency in size estimates of sections of EAS viewed 
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simultaneously at different emission angles by the two eyes. 

The size versus depth distributions are fitted with a Gaussian 
form and the X max and energy of the event determined(l) . The Gaussian 
form fits most showers well . Figure 1 shows a typical shower profile. 

3. Systematics of 
Shower Maximum Di~stri- 
butions The absol ute 
depth of X max is 
the parameter most 
sensitive to 
systematic errors of 
any that we measure. 
Symmetric random 
errors in the zenith 
angle, e z , result in 
non-symmetric errors 
in X max . Symmetric 
errors in R p also 
yield non-symmetric 
errors in X max 
because of the 
exponential atmo- 
spheric density 
distribution . 

There is also an 
Intrinsic correlation 
between X max and the shower energy. Showers whose reconstruction 
err to smaller R p and larger e z will have systematically deeper 
X max and small er estimated energy while errors that lead to larger 
R p and smaller e z yield smaller X max and larger energy estimates. 

This is a direct consequence of random errors and the exponential 
nature of the atmosphere. Any additional systematic bias in e z or 
Rp will shift the X max distributions accordingly. 

To reduce these effects to a minimum, we consider events 
that are very well reconstructed, with R p >2.0 km, projected track length 
>50°, SR p /R p <.l, Se z <10° and relative uncertainties in Gaussian 
width, Xn, ax , and energy of <.4. These cuts also have the effect of 
reducing the Cerenkov subtractions to a level where 50% variation in 
the Cerenkov light model parameters do not significantly affect the 
X max distributions. Since any residual reconstruction bias will 
affect data near the tails of the energy distribution, we cut on .7<E<3 
EeV, around the maximum of our energy acceptance. We believe residual 
systematic effects in this data sample will produce less than a ± 50 
gm/cm 2 shift in the average X max . We note that the width of the 
X m ax distribution is much less sensitive to systematic errors . We 
estimate the systematic error in the width to be ± 10 gm/cm 2 . The 
resultant distribution in X max is shown in Figure 2. The average 
Xmax for this sample is 730±60 gm/cm 2 while the width (standard 
deviation) is 120±40 gm/cm 2 . 



Figure 1. Typical Longitudinal Shower 
Profile 
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4. Monte Carlo 
Generation 
Proton and 

iron induced showers 
are generated in a 
Monte Carlo program 
with the following 
input. Protons 
are assumed to obey 
Hill as pure scaling 
with the only scale 
violating effects 
being rising cross- 
sections as 
parameterized by 
Gaisser and 
Yodh(2). The mean 
inelasticity is 
assumed to be .5. 

For the case of 
iron primaries a 
superposition model 
is assumed. The 

Monte Carlo follows hardons and electromagnetic particles down to 
l/30th of the primary energy after which parameter! zations of shower 
development are used. We use these showers to predict the number of 
photoelectrons and relative time delays observed in the detector and 
generate fake events which are then passed thru the same reconstruction 
and analysis programs as the real data. This work is in progress but 
preliminary indications are that the depth of maximum distribution for 
heavies is substantial ly narrower than that for protons. It is 
important to note that comparing the observed X max distributions to 
theoretical predictions without taking into account the details of 
detector response can lead to misleading conclusions. 

5. Conclusions. Although we are not yet ready to quote quantitative 
comparisions between the data and Monte Carlo simulations, preliminary 
evidence based on comparisons of the widths of the X m ^ x distribution 

to Monte Carlo is in favor of a substantial contribution from light 
primaries to the cosmic ray composition at these energies. Detailed 
comparisons will be forthcoming. 
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Figure 2 . Distributions of X max in 
gm/cm 2 . 




